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^ . Abstract 
(N , 

' We study the asymptotics of a Markovian system of > 3 particles in [0, 1]*^ in 

^ ^ . which, at each step in discrete time, the particle farthest from the current centre 

of mass is removed and replaced by an independent U[0, 1]'^ random particle. We 
show that the limiting configuration contains — 1 coincident particles at a random 
location ^tv € [0, 1]"^. A key tool in the analysis is a Lyapunov function based on 
the squared radius of gyration (sum of squared distances) of the points. For d = 1 
we give additional results on the distribution of the limit showing, among other 
things, that it gives positive probability to any nonempty interval subset of [0, 1], 



^ , and giving a reasonably explicit description in the smallest nontrivial case, A^ = 3. 



Keywords: Keynesian beauty contest; radius of gyration; rank-driven process; sum of 
squared distances. 
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(Secondary) 



O 1 Introduction, model, and results 

m _ 

In a Keynesian beauty contest, N players each guess a number, the winner being the 
player whose guess is closest to the mean of all the guesses; the name marks Keynes's 
^ \ discussion of "those newspaper competitions in which the competitors have to pick out the 
six prettiest faces from a hundred photographs, the prize being awarded to the competitor 
whose choice most nearly corresponds to the average preferences of the competitors as 
a whole" [71 Ch. 12, §V]. Moulin (TUl p. 72] formalized a version of the game played on 
a real interval, the "p-beauty contest" , in which the target is p (p > 0) times the mean 
value. See e.g. [2] and references therein for some recent work on game-theoretic aspects 
of such "contests" in economics. 

In this paper we study a stochastic process based on an iterated version of the game, 
in which players randomly choose a value in [0, 1], and at each step the worst performer 
(that is, the player whose guess is farthest from the mean) is replaced by a new player; 



*Department of Mathematics and Statistics, University of Strathclydc, 26 Richmond Street, Glasgow 
Gl IXH, UK. 

^Centre for Mathematical Sciences, Lund University, Box 118 SE-22100; Lund, Sweden, and Depart- 
ment of Matliematics, University of Bristol, University Walk, Bristol BS8 ITW, UK. 

■'■Department of Mathematical Sciences, Durham University, South Road, Durham DHl 3LE, UK. 



1 



each player's guess is fixed as soon as they enter the game, so a single new random value 
enters the system at each step. Analysis of this model was posed as an open problem 
in |31 p. 390]. The natural setting for our techniques is in fact a generalization in which 
the values live in [0, l]*^ and the target is the barycentre (centre of mass) of the values. 
We now formally describe the model and state our main results. 

Let d E N := {1,2,...}. We use the notation Xn = (xi, X2, . . . , x„) for a vector 
of n points Xi G M.'^. We write fini'^n) '■= ^"=i for the barycentre of Xn, and 
II ■ II for the Euclidean norm on M"^. Let ord(A'n) = X(2), . . . , denote the 

barycentric order statistics of Xi, . . . so that ||x(i) — /i„(A:'„)|| < ||x(2) — /in('^n)|| < 
■ ■ ■ < W^in) — Ain.('^n)||; any ties are broken randomly. We call X* := X(„,) the extreme 
point of Xn, a point of Xi, . . . , farthest from the barycentre. We define the core of Xn 
as A*^ := (x(i), . . . , X(n-i)), the vector of Xi, . . . , x„ with the extreme point removed. 

The Markovian model that we study is defined as follows. Fix > 3. Start with 
Xi(0), . . . , Xjv(O), distinct points in [0, 1]"', and write Xn{0) := {X^i){0), . . . , X(jv)(0)) for 
the corresponding ordered vector. One possibility is to start with a uniform random initial 
configuration, by taking Xi{0), . . . ,Xn{0) to be independent U[0, 1^ random variables; 
here and elsewhere U[0, l]*^ denotes the uniform distribution on [0, l]'^. In this uniform 
random initialization, all N points are indeed distinct with probability 1. Given X^lt), 
replace X^(t) = X(Ar)(t) by an independent U[0, 1^ random variable Ut+i, so that XN(t + 
1) = ord(X(i)(t), . . . ,X(;v-i)(t), Ut+i). 

The interesting case is when > 3: the case A^ = 1 is trivial, and the case N = 2 
is also uninteresting since at each step either point is replaced with probability 1/2 by 
a U[0, 1] variable, so that, regardless of the initial configuration, after a finite number of 
steps we will have two independent f/[0, 1] points. Our main result. Theorem II. H shows 
that for A^ > 3 all but the most extreme point of the configuration converge to a common 
limit. 

Theorem 1.1. Let d E N and N > 3. Let Xj^{0) consist of N distinct points in [0, 1]*^. 
There exists a random C,^' '■= (,n{'^n{0)) £ [0, 1]"^ such that 

X'^it) ^ i^^,^^,...,^^), and X*^{t)-Ut^O, (1.1) 

as t —7- oo. In particular, for U ~ U[0, 1^, as t ^ oo, 

^N{t) — > i^N, Cn, • • • , Cn, U). 

Remark 1.1. Under the conditions of Theorem \l.l\ despite the fact that X^{t) — Ut^O 
a.s., we will see below that X^(t) ^ Ut infinitely often a.s. 

Theorem II. II is proved in Section [2j Then, Section |3] is devoted to the one-dimensional 
case, where we obtain various additional results on the limit ^n- Finally, the Appendices, 
Sections m and El collect some results on uniform spacings and continuity of distributional 
fixed-points that we use in parts of the analysis in Section [31 

2 Proof of convergence 

Intuitively, the evolution of the process is as follows. If, on replacement of the extreme 
point, the new point is the next extreme point (measured with respect to the new centre of 
mass), then the core is unchanged. However, if the new point is not extreme, it typically 
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penetrates the core significantly, while a more extreme point is thrown out of the core, 
reducing the size of the core in some sense (we give a precise statement below). Tracking 
the evolution of the core, by following its centre of mass, one sees increasingly long periods 
of inactivity, since as the size of the core decreases changes occur less often, and moreover 
the magnitude of the changes decreases in step with the size of the core. The dynamics 
are nontrivial, but bear some resemblance to random walks with decreasing steps (see 
e.g. [3l[8] and references therein) as well as processes with reinforcement such as the Polya 
urn (see e.g. [11] for a survey). 

Our analysis will rest on a 'Lyapunov function' for the process, that is, a function 
of the configuration that possesses pertinent asymptotic properties. One may initially 
hope, for example, that the diameter of the point set X]y(t) would decrease over time, 
but this cannot be the case because the newly added point can be anywhere in [0, l]*^. 
What then about the diameter of Xj^ {t) , for which the extreme point is ignored? We will 
show later in this section that this quantity is in fact well behaved, but we have to argue 
somewhat indirectly: the diameter of A'jJ^(t) can increase (at least for big enough; see 
Remark 12.21 below). However, there is a monotone decreasing function associated with 
the process, based on the sum of squared distances of a configuration, which we will use 
as our Lyapunov function. 

For n e N and = {xi,X2, ■ ■ ■ , x„) G M*^", write 

n i—l n 

Gn{Xn) := Gn{Xi, . . . , X„) := '^""^ X] X] ~ ^j"^ ^ X] ~ Atn('^n)f ; (2.1) 

1=1 j=l i=l 

a detailed proof of the (elementary) final equality in (12. ip may be found on pp. 95-96 
of [5], for example. We remark that ^Gn is the squared radius of gyration of xi, . . . , x^: 
see e.g. [5], p. 95. Note also that calculus verifies the useful variational formula 

n 

Gn{xu...,Xn) = inf y2\\x,-yf. (2.2) 

i=l 

For n > 2, define 

Fni'^n) '■= Fn{Xi, . . . , Xn) ■= Gn-l{X^) = Gn-l{X(^i), . . . , X(„_i)). 

Lemma 2.1. Let n>2 and Xn = {xi,X2, ■ ■ ■ , x„) G M^". Then for any x G M^, 

F„(X(i), . . . ,X(„_1),X) < Fn{Xn). 

Proof. For ease of notation, we write simply (xi, . . . , x„) for (x(i), . . . , X(„)), i.e., we relabel 
so that Xj is the jth closest point to fini'^n)- Then X* = x„, X^ = (xi, . . . , x„_i), and 

n-l 

-Fold := Fn{Xn) = Gn-l{Xl, Xn-l) = ^ \\Xi - /iold||^ (2.3) 

i=l 

where /i^j^j := yU„_i(A'^). We compare Fold to F„ evaluated on the set of points obtained 
by removing Xn and replacing it with some x G M"^. 

Write y := {xi, . . . , x„_i, x}* for the new extreme point. Then 

n-l 

Fnew •= Fn{Xi, . . . ,X„,_i,x) = ^ ^ ||Xj A'-ncwII ~^ \\-^ /^ncwll 11^ ~ A'-ncwII ' (2-4) 

i=l 
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where 

'n-l 



n 



^—[Y,x, + x-y\= /i^i, + (2.5) 



Denote /incw := fin{xi, . . . , x^-i, x), so 

/^ncw = T 7- 2.6 

n — \ n — 1 

From (ESD, (El, and (ES]), we obtain 

n-l 
i=l 

For the sum on the right-hand side of (12.71) . we have that 

n— 1 n— 1 

(Iki - /inewf - Iki - /ioldf ) = XI (2^* ■ (^old - /^ncw) + IIKcwf " ll/^oldf ) 



i=l i=l 



= (n - 1) (2/i;id ■ (/i^id - /i^ew) + ll/^ncwf - ll/^oldf) 
= (n - 1) (ll^df - 2(/X:,M ■ /incw) + ll/i'newir) ■ 

Simphfying this last expression and substituting back into (12. 7p gives Fnew — -^oid = 
- l)||/ioid - /incwll^ + 11^; - Kcwf - \\y - /incwll^- Thus, usiug ^ and then 

Fnew - Fold = + ||xf - ||yf - 2Kew " " y) 

n — \ 



n — 1 \n — 1 n — 1 

Hence we conclude that 

Kew - Fold = 7 - \\yf - 2/inew " {x - y)) 

n — 1 ^ ' 



n 



|a; - /incwf - ||y - /incwD < 0, (2.^ 



n — 1 

since y is, by definition, the farthest point from /inew CH 

Consider F{t) := Fp^{Xiq{t)). Lemma [2.11 has the following immediate consequence. 
Corollary 2.1. Let N>2. Then F{t + 1) < F{t). 



Corollarv 12 . 1 1 shows that our Lyapunov function F{t) is nonincreasing; later we show 
that F{t) — > a.s. (see Lemma [2.41 below) . First, we need to relate F{t) to the diameter 
of the point set X'^{t). For n > 2 and G M , write 

Dn{,~Cli ■ ■ ■ ■, Xn) '.= max 

l<i,j<n 

Lemma 2.2. Let n >2 and Xi, . . . ,Xn G M°'. Then 

^Dn{Xi,. . . ,Xnf < Gn{Xi,...,Xn) < ^ (^^ " l)Dn{Xi, . . . , Xn)"^ ■ 
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Remark 2.1. The lower bound in Lemma \2.S\ is sharp, and is attained by collinear 
configurations with two diametrically opposed points Xi,Xj and all the other n — 2 points 
at the midpoint ^2{xi,Xj) = /Xn(a^i, • • • ,Xn)- The upper bound in Lemma \27^ is not, in 
general, sharp; determining the sharp upper bound is a nontrivial problem. The bound 
Gn{xi, . . . , Xn) < f (d+i) ^n{xi, . . . , 17^ is also not always sharp. Witsenhausen 
fTSf conjectured that the maximum is attained if and only if the points are distributed as 
evenly as possible among the vertices of a regular d-dimensional simplex of edge-length 
Dn{xi, . . . ,Xn); this conjecture was proved relatively recently fUm^- 

Proof of Lemma \2.2[ Fix Xi, . . . , x„ G M*^. For ease of notation, write /i = fin{xi, • • • , Xn)- 
First we prove the lower bound. For n > 2, using the second form of G„ in (12. ip . 

n 

Gn{xi, . . . , Xn) = ^ \\Xi - > \\Xi - /if + \\Xj - 

i=l 

where (xi,Xj) is a diameter, i.e., Dn{xi, . . . = \\xi — Xj\\. By the n = 2 case of (12. 2p . 

W^i - + ll^j - A^lP > 2||xi - fi2{xi,Xj)\\^ = ^\\xi - Xjf. 
This gives the lower bound. For the upper bound, from the first form of Gn in (12. ip . 

1 " 

Gn{Xi, ...,Xn)<- - l)Dn{Xi, Xnf , 

^ 1=1 

by the definition of D„, which yields the result. □ 

Let D{t) := DN-i{X'^{t)). 

Remark 2.2. By Lemma\M (or /fO)). G2{X^{t)) = ^D2{X^{t)f, so when N = 3, 
Lemma \2. 1\ implies that D{t + 1) < D{t) a.s. as well. If d = 1, it can be shown that D(t) 
is nonincreasing also when N = A. In general, however, D(t) can increase. 

Let J^t := (t(A7v(0), X]\f{l), . . . , Xfif{t)), the a-algebra generated by the process up to 
time t. Let B{x; r) denote the closed Euclidean (i-ball with centre x G M"^ and radius 
r > 0. Define the events 

A+i := {Ut+i G fi(/ijv-i(A';(t)); 3D(t))}, A[_,, := {f/^+i G B{fiM-i{Km 

Lemma 2.3. There is an absolute constant 7 > for which, for all N > 3 and all t, 

^Ui ^ {Fit + 1) - Fit) < -lN-'F{t)} C {F{t + 1) - F{t) < 0} C A^+i. (2.9) 

Moreover, there exist constants c > and C < 00, depending only on d, for which, for 
all N > 3 and all t, a.s., 

P [F{t + 1) - F{t) < --fN~^F{t) I Tt] > PK+i I J^t] > cA^"'^/^(F(t))^/2; (2.10) 
P [F{t + 1) - F{t) <0\Tt]< P[A+i I J't] < C{F{t))'"\ (2.11) 
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Proof. For simplicity we write Xi,...,Xn~i instead of . . . , X(Ar_i)(t) and D 

instead of D{t) = Dn-i{Xi, . . . , Xn~i)- By definition of D, there exists some i G 
{l,...,iV-l} such that - > D/2, where fi'^y^ = fiN-i{Xi, . . . , Xn~i). 

Given J-^, the event A[_^^-^^, that the new point U := Ut+i falls in B{fi'^i^] D/4), has 
probability bounded below by O^D''-, where 9^ > depends only on d. Let /^new •= 
fiNiXi, . . . ,Xisf_i, U). Suppose that A'^j^^ occurs. Then, 

ll/incw - /i^ldll =]^\\U- I^'omW ^ ^ ^ (2-12) 

since > 3. Hence, by f l2.12p and the triangle inequality, 

„ „ / „ „ , .. D D AD , , 

\\U - /incwil < W - /ioldll + ll/^ncw " /ioldll - ^ + 12 " T2"' ^'^'^^^^ 

On the other hand, by another application of the triangle inequality and (12.121) . 

D D 5D 

llA'-new — llA'-old ~ ~ IIA'ticw ~ A'-oldll — ~2 Y2 

Then, by definition, the extreme point Y := {Xi, . . . , Xjv„i, U}* satisfies 

5D 

\\Y - fl^o^W > ll/incw - XiW > — . (2.14) 

Hence from the x = U case of (12. 8p with the bounds (12.131) and (I2.14p . we conclude that 

F(t + 1) - F(t) < JL. _ (i^y^ K^;^,) < -±-^DHiAU.), (2.15) 

for all > 3; the first inclusion in (12.90 follows (with 7 = 9/72) from (12.151) together 
with the fact that, by the second inequality in Lemma [2. 2^ > 2N~^F{t). This in turn 
implies fCT?]) . using the fact that W[A[^^ \ Tt] > 9dD'^. 

Next we consider the event At+i. Using the same notation as above, we have that 

ll/inew - t^ll > ll/^old - ^11 - ll/^new " /ioldll = (1-^7) IIKld " U\\, 



by the equality in (l^TT^ . Also, for any k e {I, . . . , N - 1}, 

ll/inew - -^fcll < ll/^old - -^fell + ll/^old - /^ncwll < D + — ||^old " ^^11' 

by (I2.12P again. Combining these estimates we obtain, for any A;g{1,...,A^ — 1}, 



(l - ^) y'oM -U\\-D> ^ll/i^i, - U\\ - D, 



1 1 /^ncw ^11 II /incw Xf^ 1 1 ^ ^1 

for A^ > 3. So in particular, \\fincw — U\\ > ||/incw — ^fc|| for all A; G {1, . . . , A^ — 1} 
provided Wfi'oid ~ ^11 > 3-D, i.e., U ^ B{fi'^^^]3D). In this case, U is the extreme point 
among U,Xi, . . . , X]sf-i, i.e., 

A'i^,C{X*^{t + l) = Ut+i}. (2.16) 

In particular, on A^^-^, F{t + 1) = -F(t), and F{t + 1) < F{t) only if At+i occurs, giving 
the final inclusion in (12.91) . Since | J-^] is bounded above by CdD'^ for a constant 

Cd < oo depending only on d, (12. lip follows from the first inequality in Lemma 12.21 □ 
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Lemma 2.4. Suppose that N > 3. Then, as t ^ oo, F{t) — a.s. and in . 

Proof. Let £ > and let a := mm{t G Z+ : F{t) < e}, where Z+ := {0, 1,2,.. .}. 
Then by fl2.10p . there exists 6 > (depending on e and A^) such that, a.s., 
P [F{t + 1) - F{t) < -5 I J^t] > 5l{t < a}. Hence, since F{t + 1) - F{t) < a.s. by 
Corollary 

E [F{t + 1) - F{t) I Tt] < -5^1{t < a}. (2.17) 

By Corollary 12. H F{t) is nonnegative and nonincreasing, and hence F{t) converges a.s. 
as t — )■ oo to some nonnegative limit -F(oo); the convergence also holds in since F{t) is 
uniformly bounded. In particular, E[F(t)] — )■ E[F(oo)]. So taking expectations in (I2.17p 
and letting t — t- oo we obtain 

limsup5¥[a > t] < 0, 

t^oo 

which implies that F[a > t] — as t — t- oo. Thus a < oo a.s., which together with 
the monotonicity of F(t) (Corollarv l2.1l) implies that F(t) < e for all t sufficiently large. 
Since e > was arbitrary, the result follows. □ 

Recall the definition of At and A[ from before Lemma [2T3l Define [J-'t) stopping times 
To := and, for n G N, r„ := min{t > r„_i : At occurs}. Then F(t) < F(t — 1) can 
only occur if t = r„ for some n. Since P[y4t+i | J-'j] is bounded below by a constant times 
D(tY, it is not hard to see that, provided -D(O) > 0, At occurs infinitely often, a.s., so 
that r„ < oo for all n. 

Lemma 2.5. Let N > 3. There exists a > such that, a.s., D{Tn) < e~°" for all n 
sufficiently large. 

Proof. We have from fl2.15p and the second inequality in Lemma 12.21 that 

F{t^) - F{r,, - 1) < -5F(r„ - 1)1«), 

for some 6 > 0. Note also that, by definition of the stopping times r„,, F(r„, — 1) = F(rn-i). 
Hence, 

P[F(r„) - F(r„_i) < -5F(r„_i) | J-.„_J > P« | J-,„_J > S, 
taking S > small enough, since, using the fact that 1{At-^) = 1 a.s., 

P[< I J-.„_J = E [P[< I TrMArJ I -F.„_J = E [P[< I I , 

where by definition of At and A't, F[A'^^ \ Ar„] is uniformly positive. Since F{t+1)—F{t) < 
a.s. (by Corollary 12.11) it follows that 

E [F(r„) - F(r„„i) | < -5^F{t^^{). 

Taking expectations, we obtain E[F(r„)] < (1 — (5^)E[F(r„_i)], which implies that 
E[F(r„)] = 0{e~^^), for some c > depending on 5. Then by Markov's inequality, 
AF{Tn) > e-"-^/"^] = (^(e-'^"/^), which implies that F(rn) = O^e'^''/^), a.s., by the Borel- 
Cantelli lemma. Then the first inequality in Lemma [2.21 gives the result. □ 

Remark 2.3. The proof of Lemma \2. 51 shows that F[A'^^ \ Jv„_J is uniformly positive, 
so Levy's extension of the Borel-Cantelli lemma, with the fact that r„ < oo a.s. for all n, 
shows that A[ occurs for infinitely many t, a.s. With the proof of Lemma \2.3\ this shows 
that X^{t) ^ Ut infinitely often, as claimed in Remark [7T71 
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Now we are almost ready to complete the proof of Theorem II. 1[ We state the main 
step in the remaining argument as the first part of the the next lemma, while the second 
part of the lemma we will need in Section 1331 below. For e > 0, define the stopping time 
i^e := min{t G N : F(t) < e^}; for any e > 0, z/^ < oo a.s., by Lemma 

Lemma 2.6. Let N > 3. Then there exists G [0, 1]'^ such that /ijv-i('^iv(^)) ~^ cl-S- 
and in L"^ as t —t- oo. Moreover, there exists an absolute constant C such that for any 
e > 0, and any to G N, on {u^ < to}, a.s., 



E 



max||/iAr_i(A'^(t)) - HN-li^Ni^oM I J^to 

t>to 



< Ce. 



Proof. Let := /i7v-i(A'j^(t)). Observe that for > 3, X'f^{t) and X'^{t - 1) have 
at least one point in common; choose one such point, and call it Z{t). Then /i'(^) ^ 
hull A'jJ^r(t) C hullA:'7v(^); where huUAf denotes the convex hull of the point set X. So 
\\Z{t) - ^\t)\\ < D{t). Similarly \\Z{t) - ^'{t - 1)|| < D{t - 1). By definition of 
r„, = — 1) and D{t) = D[t — 1) unless t = Tn for some n, in which case 

/i'(r„ - 1) = fi'{Tn-i) and D(r„ - 1) = ^(r^.i). Hence, 



J2 WAt) - f^'it - 1)11 = 5^ \W{r^) - f^'MW 

t>l n>l 

< E - + - ' (2.18) 

n>l 

by the triangle inequality. Then the preceding remarks imply that 

E ll/i'W - t^'it - 1)11 < J2iDiTn) + DiTn-l)) < OO, a.S., 
t>l n>l 

by Lemma 12.51 Hence there is some (random) C,n £ [0, 1]°' for which /i'(t) — )■ C,n a-s. as 
t — OO, and convergence follows by the bounded convergence theorem. 

For the final statement in the lemma we use a variation of the preceding argument. 
Let M := max{?T, G Z+ : r„ < to}. Then F{to) = F{tm) and Ii'{tm) = l^'ito), so that on 
{^e ^ we have {v^ < tm} as well. Hence (by Corollary 12.11) F{tm) < e^- A similar 
argument to that in the proof of Lemma [2.51 shows that, for m > 0, 

E[F(rM+™) I J't,] < e"^™E[F(rM) | J't,] < e\-^^ , 

on {^e ^ ^o}) where c > depends on but not on m or e. Thus by Lemma [2.2[ on 
{ve < to}, E[L)(rA/+„,)2 I < 2e^e-^"^. Also, similarly to (HJM . 

max ||/i'(t) - /i'(r,,)f < ^ ll/^'W - f^'i^ - l)|p 

t>TM ^ ' 

t>T]\I 

< ^{D{TM+rn) + D {tm . 

m>l 

Taking expectations and using the Cauchy-Schwarz inequality, we obtain, on {u^ < to} , 
E[max||/i'(t) -/i'(to)ir | J^t<] = E [max ||/i'(t) - fi'{TM)f \ Tt^ < 8eV V e''^'", 

m>l 

which is a constant times e"^. The result follows from Jensen's inequality. □ 
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Proof of Theorem \l.l[ Again let := /iAr_i(A'j^(t)). We have from Lemma [2.61 that 

— )■ a-s. Now, for any j G {1, . . . , — 1}, by the triangle inequality, 



Xu)(t) - i^w < wxu.it) - + Wit) - < D{t) + y{t) - i 



N\ 



which tends to a.s. as t — )■ oo, since D(t) — )■ a.s. by Lemma [2.51 This establishes the 
first statement in (11.11) . Moreover, by fl2.16p . X^{t + 1) 7^ Ut+i only if At+i occurs. On 
At+i, X^{t + 1) is one of the points of Xjy{t), and so in particular + 1) — — 

D{t). In addition, on A+i, we have ||f/t+i — fi'(t)\\ < 3D(t). So by the triangle inequality. 



□ 



\mt + l)-Ut+i\\<4D{t)l{A+^), 
which tends to a.s., again by Lemma [2. 5[ This gives the final part of fll.ll) . 



3 The limit distribution in one dimension 
3.1 Overview and simulations 

Throughout this section we restrict attention to d = 1. Of interest is the distribution 
of the limit in f H-iP - and its behaviour as — > 00. Simulations suggest that is 
highly dependent on the initial configuration: Figure [1] shows histogram estimates for C,n 
from repeated simulations with a deterministic initial condition. In more detail, 10® runs 
of each simulation were performed, each starting from the same initial condition; each 
run was terminated when D[t) < 0.0001 for the first time, and the value of fiN-i{X'^{t)) 
was output as an approximation to C,n (cf Theorem 11.11) . Note that, by fl2.9p . in the 
simulations one may take the new points not f/[0, 1] but uniform on a typically much 
smaller interval, which greatly increases the rate of updates to the core configuration. 
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Figure 1: Normalized histograms each based on 10® simulations, with A^ = 3 and initial 
points |, |, f {left) and N = 7 and initial points |, k G {1,...,7} (right). 

Figure [2] shows sample results obtained with an initial condition of A^ i.i.d. U[0, 1] 
random points. Now the histograms appear much simpler, although, of course, they can 
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be viewed as mixtures of complicated multimodal histograms similar to those in Figure 
[TJ In the uniform case, it is natural to ask whether C,n converge in distribution to some 
limit distribution as — )■ oo. 
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Figure 2: Normalized histograms for 10^ simulations with random i.i.d. uniform initial 
conditions, with {top row) = 3, 10 and {bottom row) N = 50, 100. 

The form of the histograms in Figure |2] might suggest a Beta distribution (this is 
one sense in which the randomized beauty contest is "reminiscent of a Polya urn" [H 
p. 390]). An ad- hoc Kolmogorov-Smirnov analysis (see Table 13. ip suggests that the 
distributions are indeed 'close' to Beta distributions, but different enough for the match 
to be unconvincing. Simulations for large A^ are computationally intensive. We remark 
that it is not unusual for Beta or 'approximate Beta' distributions to appear as limits 
of schemes that proceed via iterated procedures on intervals: see for instance [6] and 
references therein. 

In the rest of this section we study C,n and its distribution. Our results on the limit 
distribution, in particular, leave several interesting open problems, including a precise 
description of the phenomena displayed by the simulations reported above. In Section IX^ 
we give an alternative (one might say 'phenomenological') characterization of the limit C,n, 
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N 


/3 




3 


1.256 


0.0010 


10 


1.392 


0.0016 


50 


1.509 


0.0018 


100 


1.539 


0.0019 



Table 1: is the Kolmogorov-Smirnov distance between a Beta(/3, /3) distribution and 
the empirical distribution from the samples of size 10® plotted in Figure O minimized 
over /3 in each case. 

and contrast this with an appropriate rank-driven process in the sense of [1]. In Section 
13.31 we show that the distribution of is fully supported on (0, 1) and assigns positive 
probability to any proper interval, using a construction permitting transformations of 
configurations. Finally, Section [33] is devoted to the case = 3, for which some explicit 
computations for the distribution of (in particular, its moments) are carried out. 

3.2 A characterization of the hmit 

Let 

Mt) ■.= j#{se{l,2,...,t}:X*M< f^N{XN{s))} , 

the proportion of times up to time t for which the extreme point was the leftmost point (as 
opposed to the rightmost). The next result shows that 7iN{t) converges to the (random) 
limit given by Theorem ll.lt we give the proof after some additional remarks. 

Proposition 3.1. Let d = 1 and N >3. Then limt^ooT^N{t) = C,n o-s- 

It is instructive to contrast this behaviour with a suitable rank- driven process (cf [1|). 
Namely, fix a parameter vr G (0, 1). Take N points in [0, 1], and at each step in discrete 
time replace either the leftmost point (with probability vr) or else the rightmost point 
(probability 1 — tt), independently at each step; inserted points are independent f/[0, 1] 
variables. For this process, results of [1] show that the marginal distribution of a typical 
point converges (as t — oo and then oo) to a unit point mass at vr (cf Remark 3.2 

in m). 

This leads us to one sense in which the randomized beauty contest is, to a limited 
extent, "reminiscent of a Polya urn" [H p. 390]. Recall that a Polya urn consists of an 
increasing number of balls, each of which is either red or blue; at each step in discrete 
time, a ball is drawn uniformly at random from the urn and put back into the urn together 
with an extra ball of the same colour. The stochastic process of interest is the proportion 
of red balls, say; it converges to a random limit vr', which has a Beta distribution. The 
beauty contest can be viewed as occupying a similar relation to the rank-driven process 
described above as the Polya urn process does to the simpler model in which, at each 
step, independently, either a red ball is added to the urn (with probability vr') or else a 
blue ball is added (probability 1 — vr'). 

Proof of Proposition \3.1l Given tq, ti, T2, . . ., = hnin^oo AtAr-i('^iv (''«)) is independent 
of Ut, t ^ {tq, Ti, . . .}, since, by fl2.16p . any such Ut is replaced at time t + 1. Let 5 > 0. By 
Theorem ll.H there exists a random T < oo a.s. for which maxi<j<jv_i I^tv ~ -^(i)(^)l — ^ 
for alH > T. 
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Since jjL]^{XM{t + 1)) = ^^j^iJiN-i{X'^{t)) + j^Ut+i, we have that for t >T, using the 
triangle inequahty, for any i G {1,...,A^ — 1}, 

IM^Nit + 1)) - <e+ IM'^Nit + 1)) - ^n\ 

N - 1 1 

Hence, for t > T, 

max \X^i^{t)-fiN{XN{t + l))\<h^M-Ut+i\+2e. (3.1) 
On the other hand, for t > f^N{^N{t + 1)) > ^^{^n — ^) + jjUt-\-i, so that for 

te{l,...,N-l}, 

N - 1 

f^Ni^Nit + 1)) - Ut+i > ^^(^Jv - f/m - e). (3.2) 
Suppose that Ut+i < C,n — Ke for some K G (1, oo). Then, from fl3.ip and (13. 2p . 

|/iAr(A'Ar(t + 1)) - - maX |X(i)(t) - /iAr(A'Ar(t + 1))| 

l<i<N—l 

, 3iV-l 

>£-((N-2)K-3N+l). 

This last expression is positive provided K > ^i^, which is the case for all > 3 with 
the choice K = 8, say. Hence, with this choice of K, {Ut+i < C,n — 8e} implies that Ut+i 
is farther from /i7v+i('^7v(^ + 1)) than is any of the points left over from X^(t). Write 
Lt := {Ut < — 8e}. Then we have shown that, for t > T, the event Lt implies that 
Ut = X^{t), and, moreover, Ut < ^^{XNit)). Hence, for t >T, 



S=T s=T 

s^{ro,Ti,...} 

Given Tq,Ti,..., Ug, s ^ {tq,Ti, . . .} are independent of T and ^tv- For such an s, Us 
is uniform on Is := [0,1] \ B{^n-i{'^n{s));3D{s)), and, for s > T, D{s) < 2e so that 
Is ^ [0, max{^N — Se, 0}] U [minj^Ar + Se, 1}, 1]. Hence, given s ^ {tq, ti, . . .} and s >T, 

F[Ls] = F[Us <^N-Ss\Use Is] >^N- 8e. 

Hence, considering separately the cases > 9e and < Qe, the strong law of large 
numbers implies that 



^ J2 HLs) >^N- 9e, 



t 

s^{ro,ri,...} 

for all t sufficiently large; here we have used the fact that t — T — > oo a.s. as t — t- oo and 
i^{n G Z+ : r„ < t} = o{t) a.s., which follows from (12. lip and Lemma [2.41 Since e > 
was arbitrary, it follows that liminft^oo ^~^7rjv(t) > a-S- The symmetrical argument 
considering events of the form Rt := {Ut > ^N+8e} shows that limmit^ooi^—t'^T^Nit)) > 
1 — C,N a-s., so limsupf_^oo t~%Ar(t) < a.s. Combining the two bounds gives the 
result. □ 
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3.3 The limit has full support 

In this section, we prove that is fully supported on (0, 1) in the sense that ess inf = 0, 
ess sup .^AT = 1, and C,n assigns positive probabihty to any non-null interval. Let 

m^iO) := min{||X,(0) - X,(0)|| : j G {0, 1, . . . , iV + 1}, z ^ j}, (3.3) 

where we use the conventions Xo(0) := and Xm+i{0) := 1. For p > let Sp denote 
the J-Q-event 5*^ := {mAr(O) > p} that no point of Ajv(O) is closer than distance p to any 
other point of Xi\i{0) or to either of the ends of the unit interval. 

Proposition 3.2. Let d = 1 and N > 3. Let p G (0, 1). For any non-null interval subset 
I of [0, 1], there exists 6 > (depending on N, I, and p) for which 

e I I > 61{Sp), a.s. (3.4) 

In particular, in the case where A'Ar(O) consists of N independent U[0,1] points, P[^Ar G 
/] > for any non-null interval I C [0, 1]. 

We suspect, but have not been able to prove, that has a density /at with respect 
to Lebesgue measure, i.e., ^jy is absolutely continuous in the sense that for every e > 
there exists 6 > such that P[^Ar & A] < e for every A with Lebesgue measure less than 
S. Were this so, then Proposition 13.21 would show that we may take fN{x) > for all 
X G (0, 1). Note that P[^Ar G A \ Xn{0)] may be if Xn{0) contains non-distinct points: 
e.g. if > 3 and X^^O) = {x,x, . . . ,x,y), then X^{t) = {x,x, . . . ,x) for all t. 

For a G [0, 1], e > 0, and t G N, define the event 

N 
i=l 

The main new ingredient needed to obtain Proposition 13.21 is the following result. 

Lemma 3.1. Let N > 3. For any p G (0, 1) and 5 > there exist t^ ^ N and 6q > 
(depending on N, p, and e) for which, for all a G [0, 1], 

P[E„,,(to) I > 5ol(5p), a.s. 

Proof. Fix a G [0, 1]. Let p G (0, 1) and e > 0. It suffices to suppose that e G (0, p), since 
Ea,eif) ^ Ea^^i{t) for s' > E. Supposc that Sp occurs, so that mAr(O) > p with mj^ifS) 
defined at (13. 3p . For ease of notation we list the points of Xj^iQ) in increasing order as 
< Xi < X2 < ■ ■ ■ < < 1. Let M = [N/2\ . 

Let V = e/N"^. The following argument shows how one can arrive at a configuration 
at a finite (deterministic) time to where all of Xi(to), • • • , X^ito) lie inside {a — e,a-\- e) 
with a positive (though possibly very small) probability. 

Let us call the points which are present at time old points; the points which will 
gradually replace this set will be called new points. We will first describe an event by 
which all the old points are removed and replaced by new points arranged approximately 
equidistantly in the interval [Xm, Xm+i], and then we will describe an event by which 
such a configuration can migrate to the target interval. 

Step 1. Starting from time 0, iterate the following procedure until a new point 
becomes an extreme point. The construction is such that at each step, the extreme point 
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is one of the old points, either at the extreme left or right of the configuration. At each 
step, the extreme old point is removed and replaced by a new U[0, 1] point to form the 
configuration at the next time unit. We describe an event of positive probability by 
requiring the successive new arrivals to fall in particular intervals, as follows. The first 
old point removed from the right is replaced by a new point in {Xm + i^,Xm + 1^ + 6), 
where 6 G (0, u) will be specified later. Subsequently, the ith point {i > 2) removed from 
the right is replaced by a new point in {Xm + 1^, Xm + + 5). We call this subset of new 
points the accumulation on the left. On the other hand, the ith extreme point removed 
from the left {i G N) is replaced by a new point in {Xm+i — ii'.Xm+i — iv + 5). This 
second subset of new points will be called the accumulation on the right. 

During the first M steps of this procedure, the new points are necessarily internal 
points of the configuration and so are never removed. Therefore, there will be a time ti G 
[M, A^] at which, for the first time, one of the new points becomes either the leftmost or 
rightmost point of X]s[{ti)\ suppose that it is the rightmost, since the argument in the other 
case is analogous. If at time ti the accumulation on the right is non-empty, we continue 
to perform the procedure described in Step 1, but now allowing ourselves to remove new 
points from the accumulation on the right. So we continue putting extra points on the 
accumulation on the left whenever the rightmost point is removed, and similarly putting 
extra points to the accumulation on the right whenever the leftmost point is removed, 
as described for Step 1. Eventually we will have either (a) a configuration where all the 
new points of the left or the right accumulation are completely removed, and there are 
still some of the old points left, or (b) a configuration where all old points are removed. 
The next step we describe separately for these two possibilities. 

Step 2(a). Without loss of generality, suppose that the accumulation on the right is 
empty, so the configuration consists of k points of the left accumulation and N — k old 
points remaining to the left of Xm (including Xm itself). Note that Step 1 produces at 
least M new points, so M < k < N — 1, since by assumption we have at least one old 
point remaining. Let us now denote the points of the configuration xi < X2 < ■ ■ ■ < xn 
so that XN-k = Xm, and by the construction in Step 1, XN-k+i £ {XM + ii^,XM + ii^ + S) 
for i = 1,2, . . . , k. Provided that k < N — 2, so that there are at least 2 old points, we 
will show that Xi is necessarily the extreme point of the configuration. Indeed, writing 
= l^'N^Xi, . . . ,xn), using the fact that xjq-k+i > x^^k + iu ioi 1 < i < k and xn < 
XN-k + kv + 5, we have 

xi + xn Xl^ VxN-k + kxN-k + \vk{k + l) xi+ XN^k + i^k + 5 



The old points all have separation at least p, so ioi 1 < i < N — k, xi > Xi + {i — \)p, 
and hence 



2xi + ■ ■ ■ + 2xN-k + {2k - N)xN-k > Nxi + p{N -k~l){N-k)+ p{2k -N){N -k-l). 
It follows, after simplification, that 



2 




5N). 



p- 



xi + Xn 



2 



>7^Ak{N-k-l){p-u)-6N) 
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provided 1 < A; < iV — 2. By choice of z/, we have z/ < p/9 and it follows that the last 
displayed expression is positive provided 6 is small enough compared to p (5 < p/4, say). 
Hence — /i| > \xn — fJ'\- Thus next we remove xi. We replace it similarly to the 
procedure in Step 1, but now building up the accumulation on the left. We can thus 
iterate this step, removing old points from the left and building up the accumulation on 
the left, while keeping the accumulation on the right empty, until we get just one old 
point remaining (i.e. until k = N — 1); this last old point will be Xm- At this stage, 
after a finite number of steps, we end up with a configuration where the set of points 
Xi < X2 < ■ ■ ■ < xn satisfies Xt G [Xm + {i — ^)^, + ~ + 6], i = 1,2, . . . , N. 

Step 2(b). Suppose that the configuration is such that all old points have been 
removed but both left and right accumulations are non-empty. Repeating the procedure 
of Step 1, replacing rightmost points by building the left accumulation and leftmost 
points by building the right accumulation, we will also, in a finite number of steps, 
obtain a set points Xj such that Xj G [6 + (z — l)i^, b + {i — l)v + S\, i = 1,2, . . . , N, for 
some b G [0, 1]. 

Step 3. Now we will show how one can get to the situation where all points lie inside 
the interval {a — e,a + e) starting from any configuration in which 

Xi e [b + iu,b + iu + 5], i = 1,2, . . . , N - 1, (3.5) 

where b G [0, 1] and Xi < ■ ■ ■ < x^-i are the core points of the configuration (i.e., with 
the extreme point removed). We have shown in Step 1 and Step 2 how we can achieve 
such a configuration in a finite time with a positive probability. Suppose that a > b; 
the argument for the other case is entirely analogous. We describe an event of positive 
probability by which the entire configuration can be moved to the right. 

Having just removed the extreme point, we stipulate that the new point yi belong to 
(6 + Nu — 6S,b + Nv — 55), so yi > xjy^i is the new rightmost point provided S < vjl. 
Then to ensure that x\, and not y\, is the most extreme point we need 



b + V — 



a:i H V xjv^i + yx 

N 



b + u — 



The left-hand side of the last inequality is less than —26 while the right-hand side is more 
than so the inequahty is indeed satisfied provided TV > 3. 

S 6 5 S 5S 



m I 

u b + 2v ■■■ b+{N-l)v b + Nv 



Figure 3: Schematic of a configuration at the start of Step 3. The disks represent the 
points Xi, X2, ■ ■ ■ , Xn-1 and, on the extreme right, the new point yi. 

Hence at the next step Xi is removed. Our new collection of core points is X2 < • ■ ■ < 
xn-1 < yi- We stipulate that the next new point y2 arrive in {b + {N + l)iy — 186, b + 
{N + l)z/ — 175). So again, for 6 small enough (5 < u/lS suffices), y2 > yi and the newly 
added point (yi) becomes the rightmost point in the configuration. Again, to ensure that 
the leftmost point (^2) is now the extreme one, we require 



X2 + 1/2 



iV + 3 



< 



x2-\ h xn^i + yi + y2 

N 



b + u 



N + 3' 
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The left-hand side of the last inequality is less than —85, while the right-hand side is 
more than —246/N, and so the displayed inequality is true provided > 3. 

We will repeat this process until we remove the rightmost core point present at the 
start of Step 3, namely xn-i, located in [b + {N - l)i^,b + {N - l)u + 5]. We will 
demonstrate how we can do this, in succession removing points from the left of the 
configuration and at each step replacing them by points on the right with careful choice 
of locations for the new points. We consecutively put new points yk at locations in 
intervals 

Ak:={b+{N -l + k)u-2- 3^=5, b + {N - 1 + k)u - 2 ■ 3^6 + 6), 

for k = 1,2, . . . , N — 1. We have just shown that for k = 1,2 this procedure will maintain 
the leftmost point (x^) as the extreme one. Let us show that this is true for all 1 < < 
A^ — 1, by an inductive argument. Indeed, suppose that the original points Xi,X2, ■ ■ ■ , Xk~i 
have been removed, the successive new points Uj are located in Aj, j = 1,2, . . . , k — 1, 
and that the replacement for the most recently removed point Xk~i is the new point yk- 
Place the new point yk in A^. Provided S < 4^.-^k-i_^i , Uk > Vk-i and yk is the rightmost 
point of the new configuration, while the leftmost point is G [b + vk, b + vk + 5\. Since 
> 3 we have 



Xk + Vk , , 

: < b 



u-[3''-l]6<b 



< 



2 

Xfc H h xn-1 + H h ^fc 

N 



N-1 



2(3 + 3^ + --- + 3^)(5 
^ iV 



thus ensuring that the leftmost point Xk, and not yk, is the farthest from the centre of 
mass. 

Thus, provided S < S'^u, say, we proceed to remove all the points Xk and end up 
with a new collection of points x'^, . . . , x'j^_i satisfying the property 

x[ e[b' + tiy,b' + tiy + 6'], i = l,2,...,N - 1, 

where b' = b + {N — l)v — 5' and 5' := 3^5 (> 2 ■ 3^~^5). Thus the situation is similar 
to the one in (13. 5p but with b replaced hj b' > b + v, so the whole "grid" is shifted to 
the right. Hence, provided 5 is small enough, and 5' and its subsequent analogues remain 
such that 5' < 3~^i^, we can repeat the above procedure and move points to the right 
again, etc., a finite number of times (depending on \b — a\jv) until the moment when all 
the new points are indeed in (a — e, a + e), and the probability of making all those steps 
is strictly positive. In particular, we can check that taking b < 3~'^^/'^u will suffice. 

All in all, we have performed a finite number of steps, which can be bounded above 
in terms of A^, p, and e but independently of a, and each of which required a f/[0, 1] 
variable to be placed in a small interval (of width less than 3~^z/) and so has positive 
probability, which can be bounded below in terms of A^, p, and e. So overall the desired 
transformation of the configuration has positive probability depending on A^, p, and e, 
but not on a. □ 

Proof of Proposition \3.^A Write p'{t) := p]\r^i{Xl^(t)). Let / C [0,1] be a non-null inter- 
val. We can (and do) choose a G (0, 1) and e' > such that J' := [a — e',a + e'] C J. 
Also take /" := [a — e,a + e] C /' for e = {4BC)~^N~^/'^e' , where C is the constant in 
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Lemma [2.61 and 5 > 1 is an absolute constant chosen so that e < e' /A for all > 3. Fix 
p G (0, 1). It follows from Lemma [3.11 that, for some 5o > and to £ depending on e, 

P[{/x'(to) G /"} n {D(to) < 2e} I Xi,m > 6ol{Sp). 

By Lemma O, we have that D{to) < 2e implies that F{to) < 2Ne'^ < {e'/{2BC)f, so 
that to > ^£'/(2BC)i where v. is as defined just before Lemma [2.61 Applying Lemma [2.61 
with this choice of to and with the e there equal to e'/{2BC), we obtain, by Markov's 
inequality, 



has probability at least {6o/3)l{Sp), and on this event we have I^at — a| < e + {3e'/4:) < e', 
so G /. Hence (13. 4p follows. 

For the final statement in the proposition, suppose that Ajv(O) consists of independent 
f/[0, 1] points. In this case mjv(O) defined at fl3.3l) is the minimal spacing in the induced 
partition of [0, 1] into N + 1 segments, which has the same distribution as times 
a single spacing, and in particular has density f{x) = N{N + 1)(1 — (A^ + for 
X e [0, j^] (cf Section [ID. Hence for any p E [0, ^], we have ¥[Sp] = (1 - (A^+ 
which is positive for p = say. Thus taking expectations in (13. 4 p yields the final 
statement in the proposition. □ 

3.4 Explicit calculations for = 3 

For this section we take A^ = 3, the smallest nontrivial example. In this case we can 
perform some explicit calculations to obtain information about the distribution of ^3. In 
fact, we work with a slightly modified version of the model, avoiding certain 'boundary 
effects', to ease computation. Specifically, we do not use U[0, 1] replacements but, given 
X^it), we take Ut+i to be uniform on the interval U[min X^{t) — D{t),msixX^{t) + D[t)]. 
If this interval is contained in [0, 1] for all t, this modification would have no effect on 
the value of .^3 realized (only speeding up the convergence), but the fact that now Ut+i 
might be outside [0, 1] does change the model. 

For this modified model, the argument for Theorem 11.11 follows through with minor 
changes, although we essentially reprove the conclusion of Theorem 11.11 in this case when 
we prove the following result, which gives an explicit description of the limit distribution. 

Here and subsequently '=' denotes equality in distribution. 

Proposition 3.3. Let d = 1 and N = 3 and work with the modified version of the process 
just described. Let ^3(0) consist of 3 distinct points in [0,1]. Write p := p2{X^{0)) and 
D := D2iX^iO)). There exists a random ^3 := ^siXsiO)) G R such that X.^{t) ^ (4,-^3) 
as t — 7- 00. The distribution of ^3 can be characterized via ^3 = p + DL, where L is inde- 
pendent of {p,D), L = —L, and the distribution of L is determined by the distributional 
solution to the fixed-point equation 



P max |/i'(t) - p'{to) \ < 3e'/A \ > 1/3, a.s. 



It follows that, given Xj^j^O), the event 

{/i'(to) G /"} n {D{to) < 2e} n {le^v - p'{to)\ < 3e'/4} 




with probability 1 
with probability | 
with probability | 
with probability |, 



(3.6) 
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where E[|L|^] < oo for all k, and U ~ U[0, 1]. Writing Ok := "^[L^], we have 9k = for 
odd k, and ^2 = ^'4 = |||, and 6q = In particular, 

E[e3] = E[/i], E[ef] = E[/i2] + ^E[D% and E[e|] = E[fi^] + '^-E[fiD% (3.7) 

In the case where ^3(0) contains 3 independent U[0, 1] points, E[,^|] = |, |, | /or k = 
1, 2, 3 respectively. If Xs{0) = (i, i, |), i/^en EfCf] = i, |, if, ^ /or A; = 1, 2, 3, 4. 

We give the proof of Proposition 13.31 at the end of this section. First we state one 
consequence of the fixed-point representation (13.61) . 

Proposition 3.4. L given by 113. 6\} has an absolutely continuous distribution. 
Proof. It follows from (13. 6 p that 

P [L = i] = I P [t/ (L - i) = 1] + i P [t/ (L + i) = 2] 

+ |P[f/(L-i)=0]+iP[f/ (L + i)=0]. 

The first two terms on the right-hand side of the last display are zero, by an application 
of the first part of Lemma ISTTl with X = U, Y = L± 1/2, and a = 1, 2. Also, since U > 
a.s., F[U{L ^ 1) = 0] = P[L = ±i], and, by symmetry, P[L = 1/2] = P[L = -1/2]. Thus 
we obtain 

Hence F[L = 1/2] = P[L = -1/2] = 0. 

Each term on the right-hand side of (13.61) is of the form ±^ + V{L±^) where V is an 
absolutely continuous random variable, independent of L (namely U or U/2). The final 
statement in Lemma [5.11 with the fact that F[L = ±1/2] = shows that each such term 
is absolutely continuous. Finally, Lemma [5.21 completes the proof. □ 

In principle, the characterization (13. 6p can be used to recursively determine all the 
moments E[L'^] = 9k, and the moments of .^3 may then be obtained by expanding E[,^|] = 
E,[{fi + DL)'']. However, the calculations soon become cumbersome, particularly as n and 
D are, typically, not independent: we give some distributional properties of (//, D) in the 
case of a uniform random initial condition in Section jH 

Before giving the proof of Proposition 13. 3^ we comment on some simulations. Figure 
H] shows histogram estimates for the distribution of ^3 for two initial distributions (one 
deterministic and the other uniform random), and Table [2] reports corresponding moment 
estimates, which may be compared to the theoretical values given in Proposition 13.31 In 
the uniform case, we only computed the first 3 moments analytically, namely, i, |, | as 
quoted in Proposition 13. 3j it is a curiosity that these coincide with the first 3 moments 
of the U[0, 1] distribution. 

Proof of PropositionlEM Let fi'{t) := |(X(i)(t) + X(2)(t)) and D{t) := |X(i)(t) - X(2)(t)| 
denote the mean and diameter of the core configuration, repeating our notation from 
above. 

Consider separately the events that Ut+i falls in each of the intervals [minA'3(t) — 
D{t),mmX^{t)), [mmX^{t),mmX^{t) + lD{t)), [mmX^{t) + ^^(t), max A'^(t)), 
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Figure 4: Normalized histograms based on 10® simulations of the modified = 3 model 
with fixed { — 1/2, 1/2, 100} initial condition (left) and i.i.d. U[0, 1] initial condition (right). 



k 


1 


2 


3 


4 


5 


6 


±1 core 
U[0,1] 


0.0001 
0.5000 


0.5833 
0.3333 


0.0000 
0.2500 


1.0192 
0.2029 


-0.0005 
0.1739 


3.4765 
0.1561 



Table 2: Empirical fcth moment values (to 4 decimal places) computed from the simu- 
lations in Figure m 



[max A'3(t), max A'3(t) + /^(t)], which have probabilities |, |, |, | respectively. Given 
{^'{t), D{t)), we see, for Vt+i a U[0, 1] variable, independent of {fi'(t), D{t)), 

-D{t), Vt+iD{t)) with probability \ 



ifi'it + l),D{t + l)) 



'(/i'W- 


l+Vt+1 
2 




2-Vi+i 
4 




2-Vt+i 
4 




1+Vt+i 
2 



), |\4+iL'(t)) with probability | 



Writing mk{t) = K[D(t)'' \ ^"3(0)] we obtain from the second coordinates in (13.81) 



which implies that 



(3.8) 



mk{t) = (^^(^(2 + 2-^)J D{0)'. (3.9) 

For example, mi(t) = (5/12)*D(0) and m2(t) = (1/4)*D(0)2. 

Next we show that /i'(t) converges. From (13.81) . we have that \fj,'{t + l) — fi'(t)\ < D(t), 
a.s., so to show that fi'(t) converges, it suffices to show that D(t) < 00 a.s. But this 
can be seen from essentially the same argument as Lemma 12.51 or directly from the fact 
that the sum has nonnegative terms and ^Y^^q D{t) = Xlt^o which is finite. 

Hence fi'{t) converges a.s. to some limit, ^3 say. Extending this argument a httle, we have 
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from ([3SD that \i2'{t + 1)| < + D{t), a.s., and D{t + 1) < Vt+iD{t), a.s. Hence for 

Ui, U2, ■ ■ ■ i.i.d. U[0, 1] random variables, we have D{t) < Vi - ■ ■ VtD{0) and 



t-i 



s=0 



s=l r=l 



Here Z has the so-called Dickman distribution (see e.g. [T^ §3]), which has finite moments 
of all orders. Hence E[|/i'(t)|^ | ^^(O)] is bounded independently of t, so, for any p > I, 
{fj.'{t))P is uniformly integrable, and hence limt_^oo E[(/i'(t))'= | ^"3(0)] = E[^f | ^"3(0)] for 
any k eN. 

We now want to compute the moments of ^3; by the previous argument, we can first 
work with the moments of yu'(t). Note that, from fl3.8p . 



E[(/i'(t+ 1) 
l + (-l)^' 



1 + Vt+i 



-D{t)''E 



1 A- ( — 1 )^ 9*^+1 — 1 



6 



k+ 1 

2fc+i_i 



using the fact that E[(l + 1/t+i)'^] = In particular, E[(/x'(t+l)-/i'(t))^ | X^it)] = 

for odd k, so E[/i'(t) | A:'3(0)] = /i'(0), and hence E[^3] = limt^oo E[/i'(t)] = E[^'(0)], giving 



the first statement in (13. 7p . In addition, 

E[{^^\t + l)y-{^,'{t)r\Xsit)] 

= 2fi'{t)E[fi'{t + 1) - /i'(t) I A'3(t)] + E[(/i'(t + 1) - /i'(t))2 I A'3(t)] 



16 



Ditf. 



Hence 



t-i 



E[(/i'(t))2 - (/i'(0))2 I A'3(0)] = J2E[{f,'{s + 1))2 - (^'(.))2 I X,m 



s=0 



t~l 



I^E"^2(5) 

s=0 

7 

-5:4-^(0)^ 
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as t — )■ 00, and the limit evaluates to j^D{Oy, so that E[^3] = limt_^oo E[(/i'(t))^] = 
E[(/i'(0))2] + ^E[D{Oy], giving the second statement in dXT]) . 

Write L(/i'(0), -D(O)) = ^3(^^3(0)) emphasizing the dependence on the initial configur- 
ation through /i'(0) and -D(O). Then by translation and scaling properties 



L(ai'(0),D(0)) = /i'(0) + Z}(0)L(0,l) 



(3.10) 



So we work with L := L(0, 1) (which has the initial core points at ±|). 

We will derive a fixed-point equation for L. The argument is closely related to that 
for (13.81) . Conditioning on the first replacement and using the transformation relation 
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fl3.10p . we obtain fl3.6p . From (13. 6 p we see that \L\ is stochastically dominated by 1 + 
U\L\] iterating this, similarly to the argument involving the Dickman distribution above, 
we obtain that \L\ is stochastically dominated by 1 + Z, where Z has the Dickman 
distribution, which is determined by its moments. Hence (13. 6p determines a unique 
distribution for L with E[|L|''"] < oo for all k. 



Writing (13. 6p in functional form L = "^{L), we see that by symmetry of the form of 

^, also ^(L) = - ^(-L). Hence -L = - ^(L) = ^(-L), so -L satisfies the same 

distributional fixed-point equation as does L. Hence L = — L. 
Writing 6^ '■= which we know is finite, we get 



Here 



3=0 



l + U 



2-U 



Ijk-j 



U 



k-j 



So we get 



E 



E 



2 



U 



k-j 



£=0 

i 

t=0 



1 



k-i+1 



--■ a{kj)-, 
= :b{k,j] 



9k-j{2a{k,j) + b{kJ)). 



{3.11] 



j even, j<k 

In particular, as can be seen either directly by symmetry or by an inductive argument 
using (13. lip . 9k = for odd k. For even k, one can use (13. lip recursively to find 9k, 
obtaining for example the values quoted in the proposition. 

Note that, by (IXTUj) . E[^|] = E[(/i'(0) + LL)(0))% which, on expansion, gives the final 
statement in (13. 7p . The first 3 moments in the case of the uniform initial condition follow 
from (13. 7p and Lemma l4m For the initial condition with points |, |, |, we have Z)(0) = | 
and /i'(0) = x| + (1 ~ x)| = I 4? where x is the tie-breaker random variable taking 
values or 1 each with probability i. It follows that E[/x'(0)^] = |8-^(3*^ + 5^). Then, 

using dsini), 

j=o V-^/ 



We can now compute the four moments given in the proposition. 



□ 



4 Appendix 1: Uniform spacings 

In this appendix we collect some results about uniform spacings which allow us to obtain 
distributional results about our uniform initial configurations. The basic results that we 
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build on here can be found in Section 4.2 of [13]; see the references therein for a fuller 
treatment of the theory of spacings. 

Let Ui, U2, . . . , f/„ be independent U[0, 1] points. Denote the corresponding increasing 
order statistics f/[i] < ■ ■ ■ < U[n], and define the induced spacings by Sn,i '■= U[i\ — f/[j_i], 
i = l,...,n+l, with the conventions f/[o] := and U[n+i] := 1. We collect some basic facts 
about the Sn,i- The spacings are exchangeable, and any n- vector, such as {Sn,i, • • • , Sn,n), 
has the uniform density on the simplex A„ := {(xi, . . . , Xn) G [0, 1]" : XliLi ^« — -'-}• 

We need some joint properties of up to 3 spacings. Any 3 spacings have density 
/(xi, X2, X3) = n(n — 1)(?7, — 2)(1 — xi — X2 — 2:3)""'^ on A3. We will make use of the facts 

mm{Sn,i, Sn,2} = lSn,i, {u > 1), (4.1) 
{Sn,i,mm{Sn,2,Sn,3}) = {Sn,i,^Sn,2), {n > 2)] (4.2) 
see for example Lemma 4.1 of [TB]. Finally, for any n > 1 and a > 0, /3 > 0, 



^ "'^ "'^^ r(n + 1 + a + /3) ^ ^ 

In particular E[5^ J = for keN. 

Our main application in the present paper of the results on spacings collected above 
is to obtain the following result, which we use in Section |3]H 

Lemma 4.1. Let d = 1 and N = 3. Suppose that X^i^O) consists of 3 independent U[0, 1] 
points. Then 

{fi2{xm), D2{xm)) = {{Si + iS2)c + (1 - 5i - is2){i - 0, '2S2), 

where ( is a Bernoulli random variable with F[( = 0] = P[C = 1] = 1/2. For any k G 

^t(^«)1-^" (. + l)(.!2)(. + 3) ' ^'-'^ 
E\(u (r(0W] - 4(3A:-5 + (3^+3_i)4^(.+i)) 

niMm) ] - ^k + i)ik+2)ik+3) • ^^-^^ 

So, for example, the first 3 moments of D2{X^{0)) are |, and while the first 3 
moments 0/^2(^3(0)) are \, ^, and ^. Finally, ^ii2{X'^{^)){D2{X'^{^))f\ = ^. 

Proof. The 3 points of ^^3(0) induce a partition of the interval [0, 1] into uniform spacings 
5*1, 5*2, 5*3, 5*4, enumerated left to right (for this proof we suppress the first index in the 
notation above). For ease of notation, write D := D2{X^{0)) and fi := fi2{X^{0)) for the 

duration of this proof. Then D = min{5'2,S'3} = 51/2, by fl4.ip . Moreover, min{S'2,5'3} 
is equally likely to be either 5*2 or S3. In the former case, ^ = Si + | min{S'2, S'3}, 
while in the latter case /i = 1 — 6*4 — | min{S'2, S3}. Using (14.21) . we obtain the following 
characterization of the joint distribution of fi and D. 

{fi D) = l^^'^^ ^'^2' ^^'^^ ^^^^ probability | 

' I (1 — 5*1 — jS2, 15*2) with probability |. 

In particular, E[D''] = 2~^E[S'f], which gives (Oj) by the n = 3, a = /c, ;9 = case of 
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For the moments of fi, we have from (14 .Gp that /i has the distribution of W := Si + ^82 
with probabihty 1/2 or 1 — with probabihty 1/2. So we have 



i=o ^"'^ 



-I^m; 



J' 



where li'fc := E[H^'^]. Since Wk = IE[(5'i + 182)^], we compute 



j=0 ^-^ ' ^ ' j=0 

by the n = 3,a = k — j,(3 = j, case of (14.31) . Thus we obtain 

8(1 - 4-('=+i)) 



{k + l){k + 2){k + 3)' 



It follows that 

k 



k\ 



(j + 3)\{k-j)\ 



k\ 



(j + 3)\{k-j)\ 



(-1/4^ 



j=o ^-^ ' ^ -^^ j=0 
We deduce (14.51) . after simplification, from the claim that, for any z G 

k 



k\ 



j=0 



{j + 3)\{k-j)\ 



k\ 



2z^{k + 3)\ 

Thus it remains to verify (14. 7p . To this end, note that 

k 



[z\k + 2){k + 3) + 2 - 2z{k + 3) - 2(1 - z)''+^] . (4.7) 



8iz) 



k\ 



k + 3 



k\ 



{k + 3)\ 
k\ 



fc+3 
i=0 



A; + 3 



-zY + 2" 



1 A + 3 



A; + 3 



+ z~ 



k + 3 



(1 _ ^)fc+3 + .2^^ ^ 2)(A: + 3) - ^(A; + 3) + 1 



^3(^ + 3)! _ 

which gives the claim (14. 7p . 

For the final statement in the lemma, we have from (14.61) that 



= ]f\{8, + \82){\8l)\ + ^E[(l -8,- \82){\8l)\ = ^EfS^^] = 1, 



by gD. 



□ 



We can also obtain explicit expressions for the densities of D and /i. Since D = 81/ 2, 
the density of D is foif) = 3(1 — 2r)^ for r G [0, 1/2]. In addition, fi has density given 
by 

{4r[3(l - r) - 4r] if rG [0,1/4] 

2-4r(l-r) if r G [1/4, 3/4] . (4.8) 

4(1 - r) [3r - 4(1 - r)] if r G [3/4, 1] 
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Indeed, with the representation of fi as either W ot 1 — W with probabihty 1/2 of each, 
we have F[fi < r] = ^F[W < r] + |(1 — F[W < 1 — r]). Assuming that W has a density 
fw (which indeed it has, as we will show below), we get 

Ur) = \fw{r) + lfw{l-r). (4.9) 



Using the fact that W = Si + jS2, we can use the joint distribution of (5*1, 5*2) to calculate 

rl rl-xi ^ 

F[W < r] = dxi / dx26(l — xi — X2)l{a;i + -X2 < r} 
Jo Jo 4 

pr />{4{r-2,'i))A(l-a;i) 
= / dxi / dX26(l — Xi — X2). 

Jo Jo 



After some routine calculation, we then obtain 

'l-|(l-r)3 + |(l-4r)3 if re [0,1/4] 

-( 

3^ 



F[W < r] - 
Hence W has density fw given by 
fw{r) 

Then (jM]) follows from IK9\\ . 



l-|(l-r)3 if rG [1/4, 1] 



4(1 - r)2 - 4(1 - 4r)2 if r G [0, 1 /4] 
4(1 -r)2 ifrG [1/4,1]. 



5 Appendix 2: Continuity of random variables 

In this appendix we give some results that will allow us to deduce the absolute continuity 
of certain distributions specified as solutions to fixed-point equations: specifically, we use 
these results in the proof of Proposition 13.41 The results in this section may well be 
known, but we were unable to find a reference for them in a form directly suitable for 
our application, and so we include the (short) proofs. 

Lemma 5.1. Let X and Y he independent random variables such that X has an absolutely 
continuous distribution. Then for any a ^ we have F[XY = a] =0. Morevoer, if 
F\Y = 0] = 0, then XY is an absolutely continuous random variable. 

Proof. For the moment assume that F[X < 0], ¥[X > 0], F[Y < 0], and F[Y > 0] are all 
positive. Take some < c < d. Then 

F[XY G (c, d)] = F[\og{X) + log(F) G (logc, logrf) | X > 0, F > 0]F[X > 0]F[Y > 0] 

+ P[log(-X) + log(-F) G (logc, log d) \X <0,Y < 0]F[X < 0]F[Y < 0]. 

Note that conditioning X on the event X > (or X < 0) preserves the continuity of X and 
the independence of X and Y. Then since the sum of two independent random variables 
at least one of which absolutely continuous is also absolutely continuous (see [HI Theorem 
5.9, p. 230]) we have 

P[log(X) + log(y) G (logc, logrf) I X > 0, r > 0] = / f+{x)dx, 

J log c 
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and 



P[log(-X) + log(-F) G (log c, log c/) I X < 0, F < 0] = / /„(x)dx 

J log c 

for suitable probability densities /+ and After the substitution m = e^, this yields 

■ P[X > 0]P[F > 0]/+(logu) + ¥[X < 0]P[y < 0]/_(logn) 



¥[XY G (c, d)] 



-du. 



u 



This expression is also valid if some of the probabilities for X and Y in the numerator of 
the integrand are zero. Therefore, we have, for any < c < d. 



p[xy G (c, d)] 



f{u)du, 



(5.1) 



for some function f{u) defined for m > 0. A similar argument applies to the case c < < 0; 
then (15. ip is valid for any c < d < as well, extending f{u) for strictly negative u. In 
particular, it follows that P[Xy = a] = for a 7^ 0. 

Now if F[Y = 0] = 0, then F[XY ^ 0] = 1. Then we can set /(O) = so that f l^ 
holds for all c, G M. □ 

Lemma 5.2. Suppose that a random variable L satisfies the distributional equation 

Zi with probability pi 



Zn with probability Pn, 

where n G N, J2^=iPi = ^, Pi > 0, and each Zi is an absolutely continuous random 
variable. Then L is absolutely continuous. 

Proof. Suppose Z, has a density /j. Then for any —00 < a < b < +00 we have 



P[L G (a, 6)] = 5^p,P[Z, G (a, 6)] = ^p, / f,{x)dx = / 
i=i i=i ■^'^ 



^Pifi 



X) 



1=1 



dx, 



which yields the statement of lemma. 



□ 



Acknowledgements 

Parts of this work were done at the University of Strathclyde, where the third author 
was also employed, during a couple of visits by the second author, who is grateful for the 
hospitatlity of that institution. 



References 

[1] C. Benassi and F. Malagoli, The sum of squared distances under a diameter con- 
straint, in arbitrary dimension. Arch. Math. 90 (2008) 471-480. 

[2] E. De Giorgi and S. Reimann, The a-beauty contest: Choosing numbers, thinking 
intervals. Games Econom. Behav. 64 (2008) 470-486. 



25 



[3] P. Erdos, On the smoothness properties of a family of Bernoulh convolutions, Amer. 
J. Math. 62 (1940) 180-186. 

[4] M. Grinfeld, P.A. Knight, and A.R. Wade, Rank-driven Markov processes, J. Stat. 
Phys. 146 (2012) 378-407. 

[5] B.D. Hughes, Random Walks and Random Environments; Volume 1: Random Walks., 
Clarendon Press, Oxford, 1995. 

[6] N.L. Johnson and S. Kotz, Use of moments in studies of limit distributions arising 
from iterated random subdivisions of an interval. Statist. Probab. Lett. 24 (1995) 
111-119. 

[7] J.M. Keynes, The General Theory of Employment, Interest and Money, Macmillan, 
London, 1936. 

[8] P.L. Krapivsky and S. Redner, Random walk with shrinking steps, Amer. J. Phys. 
72 (2004) 591-598. 

[9] P.A. P. Moran, An Introduction to Probability Theory, Clarendon Press, Oxford, 1968 
(paperback ed., with corrections, 2002). 

[10] H. Moulin, Game Theory for the Social Sciences, 2nd ed.. New York University 
Press, New York, 1986. 

[11] R. Pemantle, A survey of random processes with reinforcement, Probab. Surv. 4 
(2007) 1-79. 

[12] M.D. Penrose and A.R. Wade, Random minimal directed spanning trees and 
Dickman-type distributions. Adv. in Appl. Probab. 36 (2004) 691-714. 

[13] M.D. Penrose and A.R. Wade, Limit theory for the random on-line nearest-neighbor 
graph. Random Structures Algorithms 32 (2008) 125-156. 

[14] F. Pillichshammer, On the sum of squared distances in the Euclidean plane. Arch. 
Math. 74 (2000) 472-480. 

[15] H.S. Witsenhausen, On the maximum of the sum of squared distances under a dia- 
meter constraint, Amer. Math. Monthly 81 (1974) 1100-1101. 



26 



